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AN  OPERATIONAL  DECISION  MODEL  EMPLOYING 
OPERATIONAL  AND  ENVIRONMENTAL  FACTORS 


I.  Introduction  > 

The  meteorologist  engaged  in  operational  forecasting  is  constantly  striving  to  improve 
his  forecast.  Nevertheless,  the  ultimate  in  forecasting  shill,  that  of  predicting  all  future 
meteorological  events  exactly  correct,  is  beyond  him  now  and  will  “main  beyond  his  reach  in 
the  foreseeable  future.  Thus,  in  all  weather  forecasts  there  will  r-.main  a  certain  degree  of 
uncertainty.  Whether  this  uncertainty  is  a  hindrance  to  the  user  of  the  weather  information 
or  not  is  highly  dependent  upon  the  type  of  activity  being  planned  and,  consequently,  in  most 
cases  is  quite  variable.  (The  thesis  of  this  Technical  Note  is  the  demonstration  of  a  means 
which  will  provide  a  measure  of  the  degree  of  usefulness  of  meteorological  information  under 
various  situations^} 

The  measure  of  the  usefulness  and  also  the  effectiveness  of  weather  information  can  be 
determined  by  applying  well-known  rules  of  decision  theory  to  operational  situations  which  are 
affected  in  some  way  by  meteorological  phenomena.  This  approach  is  not  new  and  an  excellent 
historical  development  can  be  found  in  a  paper  by  Glahn  [1],  The  pioneering  efforts  of  Glahn 
did  not  go  unnoticed.  Cummings  [2]  pointed  out  the  value  and  need  of  such  an  approach.  Sub¬ 
sequently,  Huschke  and  Rapp  [3]  performed  a  similar  analysis  of  a  broadscale,  relatively - 
complex,  military  operation.  In  another  study  Huschke  [h]  discussed  the  use  of  meteorological 
decision  theory  in  planning  requirements  and  specifications  of  weapons  systems.  Rapp  [5] 
developed  a  simple  model  to  stress  the  usefulness  of  weather  information  in  activities  ranging 
from  the  simplest  to  the  most  complex. 

The  information  to  be  gained  from  analyses  such  as  those  mentioned  above  is  invaluable  to 
the  decision-maker.  In  the  same  respect,  the  AWS  forecaster  must  impart  es  much  pertinent 
information  as  possible  to  his  customer  and,  because  of  this,  it  is  imperative  that  procedures 
like  those  discussed  here  are  not  overlooked.  Indeed,  an  analysis  of  this  type  is  an  excel¬ 
lent  vehicle  to  further  enhance  AWS-customer  relationships  since  it  requires  strong  inter¬ 
action  between  AWS  personnel  and  user  agencies.  Active  user  participation  in  this  type  of 
analysis  is  a  prerequisite  to  obtaining  meaningful  results.  (This  requirement  has  been  re¬ 
peatedly  stressed  by  Huschke  [6]  in  establishing  guidelines  for  studies  of  this  type.) 

II.  The  Concept. 

In  actuality,  most,  if  not  all,  inputs  to  a  decision  are  uncertain;  it  is  the  degree  and 
quantification  of  this  uncertainty  that  provides  useful  information  to  the  decision-maker, 
furthermore,  measurements  of  uncertainties  must  be  in  terms  that  the  operational  decision¬ 
maker  can  easily  understand  and  use. 

In  evaluating  the  effectiveness  of  meteorological  information,  the  principal  uncertp.iuty 
must  necessarily  be  the  meteorological  information  itself  and  the  task  of  quantifying  that  in¬ 
formation  belongs  to  the  meteorological  community. 

With  this  in  mind,  the  basic  concept  is  to  model  various  operational  environments  with 
weather  conditions  or  parameters  as  tne  sole  inhibiting  factor  (in  a  "simple1'  model).  It  must 
be  assumed  at  the  outset  that  decisions  regarding  other  variables  within  the  problem  are 


1 


USAFETAC  TN  72-8 


November  1972 


correct.  Thus,  the  end  result  Is  a  true  measure  of  the  contribution  of  meteorological  informa¬ 
tion. 

A  simple  model  will  adequately  serve  to  exemplify  the  uses  and  also  the  importance  of  a 
decision  theory  approach.  Rather  than  use  an  actual  operational  situation,  it  was  considered 
more  feasible  to  contrive  a  simple,  but  somewhat  realistic,  situation  for  illustrative  purposes. 
Even  so,  the  meteorological  data  used  in  this  model  were  obtained  from  actual  climatological 
records. 

III.  The  Procedure 

The  basic  methodology  employed  in  the  ensuing  decision -model  is  that  espoused  in  a  book 
on  decision  theory  by  Chernoff  and  Moses  [7].  By  combining  conditional  probabilities,  clima¬ 
tological  probabilities,  and  operational  loss  values  for  specified  actions  in  the  manner  to  be 
described  below,  one  can  obtain  "optimum"  operational  strategies.  Conditional  probabilities 
comprise  the  "experimental  matrix".  This  matrix  contains  probabilities  that  a  particular  fore¬ 
cast  category  (F^)  for  some  meteorological  variable  will  have  been  forecast,  given  that  a 
particular  observational  category  (#^)  occurs.  These  probabilities  are  represented  in  the 
text  by  P(Fk|0i).  Climatological  probabilities,  collectively  termed  the  "climatology  vector", 
represent  the  probabilities  of  occurrence  of  the  observational  categories  based  on  past  records. 


These  probabilities  are  represented  in  the  text  by  P( 0^ ) .  The  "loss  table"  contains  opera¬ 
tional  costs  incurred  for  particular  actions  (a.)  when  specified  observational  categories 
occur.  These  values  are  represented  by  L( O^Aj).  All  of  the  above  values  are  combined  to  ob¬ 
tain  expected  (long-run  average)  losses  by  observational  category,  L^, Sq),  and,  subsequently, 
the  expected  loss  for  each  of  the  possible  strategies,  L(S^).  The  strategy  which  is  opera¬ 
tionally  feasible  and  shows  the  minimum  106s  is  then  selected  as  the  optimum  strategy.  In  this 
case,  a  strategy  is  defined  as  a  vector  (S^)  composed  of  the  actions  taken  corresponding  to 
each  forecast  category.  Thus,  S  =  (A  (F  ),  A  (F?),  .  .  .  ,  A  (F  )).  S  is  equivalent  to 
"action  order"  in  the  attached  computer  program.  For  the  above  variables,  the  subscript 


limits  are  as  follows 


symbolic  and  conceptual  representation 


:  i  =  1,  .  .  .,  m;J=  1,  .  .  .,  n;  k  =  1,  .  .  .  rj  q  =  1, 
al  representation  of  the  above  decision-model  is  shown  in  F 


,  .  .  . ,  t.  A 
Figure  1.  The 


actual  procedure  for  computing  the  L(S^)  Is  discussed  below. 

In  many  Instances  the  use  of  a  computer  would  be  necessary  for  computing  the  L(S^)  since 
many  strategies  can  result  from  having  only  a  few  forecast  categories  and  a  few  actions.  The 
total  number  of  possible  strategies  for  a  given  problem  is  equal  to  the  number  of  possible 
actions  raised  to  a  power  equaling  the  number  of  forecast  categories.  Thus,  if  there  are  4 
possible  actions  and  3  forecast  categories,  there  are  4s  or  64  strategies  that  must  be  simu¬ 
lated.  If  a  computer  is  used  to  calculate  the  losses,  large  numbers  of  strategies  can  be  ex¬ 
amined. 

In  their  study  Chernoff  and  Moses  [7]  approach  the  decision  problem  in  several  ways.  The 
technique  presented  in  this  technical  note  is  but  one  of  their  methods,  however,  it  is  felt  to 
be  an  extremely  useful  one.  Even  so,  more  investigation  into  the  field  of  decision  theory  must 
be  accomplished  before  the  best  method  can  be  determine^.  for  each  specific  type  of  problem. 
Elaboration  on  the  method  chosen  for  this  study  follows. 

A  single  procedure  is  followed  to  determine  the  expected  losse„  for  all  strategies  in  a 
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given  problem.  The  procedure  takes  into  consideration  the  cost  of  taking  specified  actions 
under  specified  observational  categories  and  the  probabilities  associated  with  the  occurrence 
of  the  observational  categories  and  the  forecast  categories.  The  first  procedural  step  to 
compute  the  expected  loss  for  a  given  strategy,  S^,  is  to  determine  the  conditional  probability 
for  each  action  with  respect  to  each  observational  category,  P(Aj  |i^).  Thi6  is  done  by  summing 
the  conditional  probabilities  for  each  occurrence  of  a  particular  A.  in  the  strategy  vector, 

S^.  Probabilities  of  occurrence  for  each  action  are  thus  obtained  for  each  observational 
category.  The  second  step  is  to  multiply  the  loss  value  for  each  action/observational- 
category  combination  by  the  respective  action-probability-of -occurrence  value  obtained  in  the 
first  step.  The  third  step  is  to  multiply  each  result  from  the  second  step  by  the  appropriate 
observational-category-probability-of -occurrence  value  (climatological  probability).  The  final 
a  :p  is  to  sum  the  results  of  the  third  step  to  obtain  the  expected  value  of  the  loss  that 
would  result  if  were  implemented.  (The  above  procedure  is  presented  in  symbolic  fashion  in 
Figure  2.)  This  procedure  is  subsequently  followed  for  each  S^,  whereupon  that  strategy  which 
is  feasible  and  has  the  minimum  expected  loss  is  selected  as  the  optimum  strategy. 


STEP  1:  Compute  P^(Aj  | )  for  all  i's  and  J’s; 


where  Pq(Aj  1^) 


k  V 


with  f  J=1  if  h  »  J 

1 

*•  J*0  if  h  ^=-J 

where  A^  is  that  action  taken  where  Is  given  for 
Strategy  Vector  S^. 

STEP  2:  Compute  L(6  ,  A  )  P  (A  |0  )  for  all  i's  and  J's 
*  J  J  ^ 

STEP  3:  Compute  L^,  S^)  for  all  i'f.  and  J's; 

where  1,(0^  Sq)  =  P(»1)  L(6i,  A^)  pq(Aj  1 0i ) 

STEP  4;  Compute  L(S  ); 

m 

where  L(S  )  =  ^  4(0^,  S^) 
i=l 
m  n 

■  L  ite,,  a  )  p  (*  !«,) 

i-i  j.i 

Figure  2 

Procedure  for  Computing  Expected  Loss  for  Strategy 
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The  above  technique  is  best  iilus ’rated  in  the  use  of  a  simple,  but  somewhat  realistic, 
operational -decision  problem.  First,  a  scenario  for  the  problem  is  presented,  followed  by  a 
listing  of  the  operational  and  environmental  data  used  for  the  problem.  A  dry-run  will  then 
be  performed  on  two  of  the  strategies  to  illustrate  the  above  procedure.  Following  this,  a 
description  of  a  computer  program  which  applies  the  above  principles  will  be  preseuted. 

IV.  Scenario  and  Dry-Run. 

The  scenario  will  consist  of  a  logistical  mission  to  resupply  Base  Aircraft  are  to 
depart  from  Base  X  with  the  intention  of  resupplying  Base  Y.  The  surface  visibility  at  Base 
X  is  considered  to  be  the  critical  factor  in  determining  whether  or  not  the  aircraft  should 
depart  for  Base  Y.  Weather  at  Base  Y  is  assumed  to  be  good  in  this  simplified  model  (in  other 
words,  weather  at  Base  Y  is  not  considered  to  be  a  significant  factor  in  the  mission  comple¬ 
tion,  e.g.,  supplies  can  be  paradropped  into  Base  Y).  The  aircraft  may  rely  on  only  06Z  and 
122  forecasts  of  visibility  for  Base  X  (made  at  00Z  and  06z,  respectively).  The  observation 
times  are  also  06z  and  12Z.  The  forecast  and  observational  categories  also  include  the  same 
visibility  intervals(a  3  miles,  <  3  miles).  It  is  risky  for  the  aircraft  to  take  off  when  the 
visibility  is  less  than  3  miles.  In  addition,  the  planes  must  depart  at  either  ObZ  or  12Z. 
There  are  three  possible  actions  that  can  be  taken  as  a  result  of  a  particular  forecast. 

These  actions  are  as  follows:  A^  =  takeoff  —  requires  2  hours  preparation  time  prior  to 
departure;  Ag  =  cancel  mission  —  i.e.,  mission  failure;  A^  =  delay  6  hours. 

Operati-  \al  losses  are  distributed  as  follows:  A  loss  of  1  unit  will  be  incurred  for 
each  hour  of  delay  (takes  into  consideration  manhours  lost  to  idle-time,  idle  planes  an* 
equipment,  etc.).  A  loss  of  2  units  will  be  incurred  for  each  hour  of  preparation  time  re¬ 
quired  (for  fuel,  maintenance,  man-hours,  etc.).  A  loss  of  7  units  will  occur  for  each 
mission  completion  (due  to  fuel  consumed  enroute,  wear  and  tear  on  equipment,  etc.).  A  lo’s 
of  5  units  will  result  from  a  wrong  decision  (e.g.,  action  A^  =  takeoff  is  taken  when  observa¬ 
tional  category  ^  =  <  3  miles  occurs).  Finally,  a  loss  of  20  units  is  incurred  for  failing 
to  complete  the  mission  (this  value  may  be  a  quantification  of  the  need  for  supplies  at  Base 
X). 

Figure  3  contains  the  loss  table,  climatology  vector,  and  the  experimental  matrix  for 
the  sample  problem.  For  this  example,  the  experimental  matrix  has  been  constructed  using  a 
persistence  forecast  (i.e.,  persistence  probability  values  were  used  to  construct  the  matrix). 
The  loss  taule  was  determined  in  a  somewhat  subjective  manner  by  the  authors,  and  assumes  that 
the  decision-maker  will  follow  the  chosen  strategy  regardless  of  the  consequences.  Remember, 
the  object  is  to  pick  that  strategy  which  will  result  in  the  lowest  long-run  average  loss. 

The  best  strategy  may  backfire  on  occasion,  but  it  should  do  so  less  often  than  tW-  other 
possible  strategies.  The  climatology  vector  and  the  experimental  matrix  arc,  or  have  been 
derived  from,  real  data.  The  loss  table  was  constructed  by  summing  the  possible  losses 
according  to  the  criteria  cited  in  the  previous  paragraph  with  respect  to  the  possible  com¬ 
binations  of  action  and  observational  category.  For  example,  if  the  observational  category 
is  <  3  miles,  and  the  action  implemented  is  a  takeoff,  considerable  cost  could  be  incurred. 
There  will  be  a  loss  of  4  units  for  2  hours  of  preparation  time.  There  will  also  be  a  5-unit 
loss  for  making  a  wrong  decision,  and  a  20-unit  loss  for  failing  to  complete  tne  mission. 
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TAKEOFF  CANCEL  OELAY 


<3  mi.  6i 


5  3  mi.  0, 
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.31 


.69 


LOSS  TABLE 


CLIMATOLOGY 

VECTOR 


<3ml.  53mi 


<3  mi.  0, 


>3  ml.  ©2 


.47 

.53 

.05 

.95 

EXPERIMENTAL 

MATRIX 


Ur,  and  Fj  volid  tor  06  Z  but  forceotf  of  00  Zj  0,  and  0,  valid  for  06  Z) 


Figure  3.  Data  for  Sample  Problems 

Thus,  the  total  loss  for  this  combination  of  action  and  observational  category  is  29  units. 

The  set  of  tables  and  data  shown  in  Figure  3  is  valid  for  the  OOZ  decision  time  (i.e.,  a 
decision  as  to  the  action  to  be  taken  is  made  &b  a  result  of  the  forecast  it  sued  at  OOZ.  An 
extension  oi  the  model  would  be  to  recycle  the  procedures  stated  above  in  cc  n Junction  with  a 
new  set  of  data  for  the  06Z  decision  time  as  would  be  required  in  the  case  of  a  "delay"  action 
from  the  first  cycle.  The  method  of  combining  costs  for  the  two  decision  periods,  in  the  case 
of  a  delay  from  the  first  period, may  seem  obvious  at  first  glance  but  fijrth<  r  consideration 
reveals  this  not  to  be  true.  Since  it  is  not  clear  to  the  authors  how  to  best  combine  the  two 
sets  of  data  t-.  obtain  viable  results,  this  Note  contains  an  exposition  of  the  procedure  and 
results  when  dealing  with  one  cycle  only. 

Given  the  above  data  and  assumptions,  a  dry  run  will  now  oe  undertaken  for  the  first  and 
eighth  strategies  in  this  example  to  illustrate  the  technique  more  clearly.  Since,  in  this 
problem,  there  are  three  actions  and  two  forecast  categories,  we  have  3*  or  9  possible  strate¬ 
gies.  By  performing  an  ordered  permutation,  it  can  be  seen  that  =  (A^A^)  and  Sg  =  (AyAg). 
According  tc  the  procedure  outlined  above,  these  strategies  can  be  represented  as  follows: 
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l(s1)  =  1(0^  s1)  +  l{92,  s1) 

=  H\)  [I»( 8j»  a1)p1  (A1|e1)  +  L(0r  a2)p1  (A2|e1)  +  L{9V  a3)p1  (a3  |e1)3  + 

P(02)  [L(02,  A1)P1  <A1  j02)  +  L(t2,  A2)Px  (A2(02)  +  L(02,  A3)P1  (A^)] 

L(Sg)  =  L(0X,  Sg)  +  Sg) 

=  P(01)  (L(0X,  A1)Pg  (A1|01)  +  L(0r,  A2)Pg  (A2|01)  +  L^,  A3)PQ  (A^)]  + 

P(02)  [L(02,  A1)Pg  (A1|02)  +  L(02,  A2)Pq  (A2|02)  +  L(0g,  A^Pg  (A3 1«2)  ] 

Referring  to  Figure  3  and  ';he  above  equations,  L(S^)  and  L(Sg)  can  be  computed  as  follows: 

L(S  )  =  .31  x  [29X1  +  20X0  +  6X0]  +  .69  X  [llXl  +  25x0  +  lixo]  =  16. $8  units 

L(Sg)  =  .31  X  [29X0  +  20X.53  +  6X  .U?]  +  ,69X[11X0  25X.95  +  11X.05)  =  20,93  units 

Each  of  the  other  strategies  is  computed  in  the  same  manner  with  the  lowest  cost,  operation¬ 
ally  -feasib?  e  strategy  being  selected  as  the  optimum. 

V.  Computer  Program. 

As  mentioned  previously,  a  computer  program  has  been  written  which  performs  the  above 
calculations.  The  program  is  written  in  FORTRAN  IV  for  use  by  personnel  not  especially  program- 
oriented.  It  was  designed  and  programmed  on  a  DEC  PDP-10  which  was  accessed  in  a  time-sharing 
mode  using  a  Teletype  33  terminsl.  The  program  operates  in  a  "conversational"  status,  in  that 
it  allows  the  user  to  "converse"  with  the  program  as  to  how  the  program  works,  wbat  data  are 
required  for  input,  what  options  are  available,  and  what  type  of  output  is  available.  A  com¬ 
plete  listing  of  the  program,  including  appropriate  documentation  is  contained  in  Attachment  1. 

A  sample  execution  run  of  the  program  is  contained  in  Attachment  2. 

VI.  Conclusion. 

This  program  could  provide  operations  or  planning  personnel  with  the  capability  of  obtain¬ 
ing  pertinent  decision-making  information  in  a  short  period  of  time,  provided  the  computer 
capability  is  readily  available.  If  the  program  were  implemented  in  conjunction  with  a 
cathode  ray  tube  (CRT)  display  terminal,  even  faster  response  and  greater  effectiveness  of  the 
decision-making  tool  could  be  realized.  The  present  program  is  capable  of  handling  up  to 
3, strategies  (56),  but  could  be  easily  expanded  to  include  more  meteorological  variables 
and  more  diversified  actions;  thus,  providing  the  user  with  greater  flexibility.  However,  the 
run  times  and  other  costs  for  enhancements  of  this  type  are  definitely  increased  and  a  cost/ 
benefit  analysis  would  have  to  be  performed  to  determine  the  feasibility  of  the  more  complica¬ 
ted  models.  For  example,  a  problem  that  contains  10  possible  actions  and  b  forecast  categories 
would  require  that  l,000,0vo  strategies  be  simulated,  which  could  run  a  considerable  length  of 
time,  depending  on  the  type  of  computer  h&rdvare  used. 

The  concept  of  combining  the  operaticial  and  meteorological  information  is  the  important 
factor  considered  in  this  study.  The  specific  dec is ion -theory  approach  presented  here  may  or 
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may  not  be  the  best  way  to  go.  It  is,  however,  a  step  toward  the  ultimate  goal  of  minimizing 
loss  of  resources  over  the  long-run.  There  is  a  definite  need  for  more  study  in  this  area  of 
improving  operational  decision-making  through  the  interaction  of  environmental  and  operational 
factors. 
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2.  Sample  Run 
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0  0,00000  OOO'OOOOOOOOOOOOOOOOOO  O  OOOOO  00.00  o' o 
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•C  ***###*»*#####*#*#, **1**#**.#;**#**#  MlNLOS  a**#####*#***#**##***** #***  **  #» 

THIS  PROGRAM  WILL  COMPUTE  LONG-PUN  MINIMUM  AVERAGE  LOSSES  FOR 
SPECIFIC:  OPERATIONAL  STRATEGIES'#  BASFQ  ON  OPERATIONAL  LOSS  VALUES  d 
WEATHER  OCCURRENCE  RROEAP IL  I  TT.fcS.o  THE  INPUT/CUTPUT  DOCUMENTATION  IS 
PRESENTED  IN  THE  PROGRAM  EXECUTION  PHASE  ANn  IS  GEARED  EITHER-' TO  THE 
EXPERIENCED  OR  NON-FYPERIENCFD  '1S,ER  -OF  THIS  PROGRAM  #  AS  APPROPRIATE  o 
LOGIC  DOCUMPNTATI ON  FOLLOWS*  FOP  ASSOCIATES  INFORMATION  REFER  TO 
FILES  IN  THE  SPECIAL  PROJECTS  SECTION* 

°ETtR  HALL/JUNE  ?°#197 2 


«Mc3*  M  *  )*M  #  ¥#*##«»#:»*»  a#**#**##*  TsEF.I  N  I  T  I  ONS  a***#****##*#**##*#*** 

ACNLOS  -  TEMPORARY  STORAGE  FOR  THE  LOSSFS  BY  ACTION. 

AGAIN  -  STORES  A  »Y  =  S»  OR  ‘NO'#  FOP  RECYCLE.  OF  E>TCJTION 
APRIOR  -  A  FR I  OR  I  (CLIMATOLOGICAL)  PROBABILITY  *  RR  A  Y 
EXPNCE  -  STORES  A  ’YES*  OR  •  N C TO  SEE:  WHETHER  EX PFR-LENCE,0. 

EXPYAE  -  EXPERIMENTAL. TABLE  ARRAY  FOR  CONDITIONAL  PROD  AGILITIES 
1  -  DO-LOOP  COUNTER  AND  INDEX 

IFOPMT  -  ALPHANUMERIC  ARRAY  FOR  STORING  VARIABLE  FOPMATS 
IC5ACN  -  ARRAY  THAT  STORES  ACTION  ORDER  FOR  EACH  STRATc6Y 
J  -  DO-LOOP  COUNTER  AND  INDEX 

KFOPMT  -  ARRAY  .FOR  STORING  IFO.RMT  (VARIABLE  FORMAT)  3ETMG  USED 
LOSYAB  -  LOSS  TABLE  ARRAY 

MINSRT  -  STRATEGY  NUMBER  ARRAY  WHEN  ASK  FOR  SORTED  OUT3UT 

N  -  00-L0.0P-  XOL'NTFR  AND  INDEX#  USED  IN  SORT  ROUTINE 

NACTNS  -  NUMBER  OF  ACTIONS  IN  PROBLEM 
NDRACN  -  ACTION  NUMPEP  INDEX 

NEROUT  -  NUMBER  OF  SORTED  STRATEGIES  10  BE  SOFTED  AND  LISTED 
KN  -  DO-LOOP  COUNTER  AND  INDEX#  USED  IN  SORT  ROUTINE 

NO  -  OO-LOOP  COUNTER  /  NT  iNDEX  FOR  NUMBER  OF  FORECAST 

CATEGORIES 

NOSKWO  "  INDEX  THAT  REVERSES  ORDER  OF  IC3A.CN  ARP  AY 
NOBS  -  NUMBER  OF  F-ORDCAS-T  CATEGORIES  IN  PROBLFM 
NS  -  DO-LOO?  COUNTER  AND  INDEX  FOR  NUMBER  OF  OBSERVATION 

CATEGOR  IF? 

NSTATE  -  NUMBER  OF  OBSERVATION  CATEGORIES  IN  'PR09LFM 
NSTRAT  -  NUM3ER  )F  STRATEGIES  POSSIBLE.  IN  PROELFM  • 

NWT  -  00-L0.0?  COUNTER  ANT  INDEX  USED  TO  STORE  A. FORMAT  STATEMENT 
IN  < FORM T 

C  SORTED  -  STORES  A  'YES«  OR  «NC  FOR  SORTING  STRATEGIES  BY  LOSS  ' 

C  STAL0S  -  ARRAY  FOR  CUMULATIVE  LOSS  IN  OBSERVATION  CATEGORY 

C  STRL03  -  ARRAY  FOR  UNSORTED  STRATEGY  LOSSES 

C  STRSRT  -  ARRAY  FOR  SORTFO  'STRATEGY  LOSSES 

C  YES  -  ALPHANUMERIC  ‘YES’  FOR  TESTING  ANSWERS  TO  DUFSTICNS 

C 

C 

C.  I*##*#***#***#*#*####**##**  FROGRAH  LISTING  a**#*'##**  #***  **#**>••****  *x  * 


ATCH  1 


rr\i  t  1*1 


;  <UEAf£TAC  >minlos*_P4»-A 


*?RI 


11-MIG-7?  5't  1 8  AM 


C" 

RFAL  LGbTAB' 

D.lhE'NS  XU M  LOST  AB  (  r  #E  V#EX  PTA6  ( *  #  5  >  *  I  OS  ACN.C" )  #  ®T  AL  OS  <  5  > 
r»IMt?NSiGN  iPB4;OF:X,:)#STRLOF(31?50  iSTfiSRTC^O)  jM!NSPT(50) 

’DIMENSION!  IF  OF  Hi  (  r#4)  jKFO  BM:TCr.) 

DATA  IFCRNT/^HCa’Xj^HIJi^j^HXiPTljTH^V-iF  # 

1  J>feHI3-jl2i,:H.<j3Uj,5Hi6XiFji;H6#?)  i 

2  KH(2?.X##RHI3i#il#‘5HX#«U#9H#6Y.#F#5H6*?)  # 

3  SHC25y>,*Hl3#li#s.HX#*UiEH#fcXjF#SH6oR>  / 

DATA  YES/THYFS/ 

c 

c  »»***♦»  REMARKS  jEXPLANATIONi  AND  QUESTIONS  FOR  INPUT* 

C 

TY?E  iO'Ol 

i  001  FORMAT  <»  O'  »20X#  '**#»*#  BAYES  STRATEGIES  h «»<»*♦* h  »  /'OARE  Y 0 •  > 

1  '  EXPERIENCED?  (YES  OR  N0>'//> 

ACCEPT  1002#  EXPUCS 
100  2'  FORMAT  !(AP.) 

IF  <fiXPNCF#EQoYES)  GO  T.O  *0 
T-Ype  10  0  0 

1000  FORMAT  ( « O' # 5X # ' TH  IS  PROGRAM  WILL  COMPUTE  LONG-RUN  MINIMUM  AVER*  # 

1  'AGE  LOSSES  FOR*/*  SPECIFIC  OPERATIONAL  STRATEGIES  ("BAYFS  STR'# 
,?  'ATEG  IES ")  RASED  ON  OPERA-*/*  IIONAL  LOSS  VALUES  AND  WEATHER*# 

2  '  OCC  URRFNCF  PRO?  ACTLIT  IFSo  THE  IN»UT  •  /•  PAR AMCT£RS  ARE.  AS'# 

<i  '  FOLLOWS:  (1)  VHE  NUMDES  OF  ACTIONS  AVAILABLE#'/*  (?)  THE  NUN'# 
P  *8f.R  OF  03SER  VAT!  ON  CATEGORIES  USED#  C3>  THE  NUMBER  OF'/'  FORE'# 

6  'CAST  CATEGORIES  USED#  («>  THE  OPERATIONAL  COST  (LOSS)  FOR*# 

7  '  EACH'/'  COMBINATION  OF  ACT  I ONG  ANO  OeSF R V AT  1 0 N  CA TE GO RI E$ #  '  # 

8  '  <P)  THF  CLIMATOLOG-'/ '  ICAL  PROBABILITY  OF  OCCURRENCE  OF'# 

9  '  EACH  OBSERVATION  CATEGORY#  AND'/  »  (6)  THE  COMCIH!  O.MAL  PROD  A5'  # 
A  '3ILTTY  FOR  EACH  COM®  XNATION  OF  OBSERVATION'/'  AND-  FORECAST  CA',# 
1  'TEGORIF?  U*E*#  T.Hr  PROBABILITY  THAT  A  PARTCULAR  FOR£-'V  CA'# 
2 ' ST  WAS  ISSUED. GIVEN  THAT  A  PARTICULAR  OBSERVATION  CATEGORY'/ 

3  »  OCC  URRFO) o ' //»  THF  INPUT  FORMATS  ARE  AS  FOLLOWS:  (1)#  (2»# 

«  ')#  AND  (3)  SHOULD  9Fi/«  T?Y®FO  AS  ONE  .DIGIT  INTEGERS;  («)#  (5)'# 

5  '#  AND  ( £ )  SHOULD  GK  TYPED  AG'/'  THREE  DIGIT  REAL  NUMBERS  ( INC.'  # 

6  'LUDING  DECIMAL  POINT  INSERTED  AS'/'  APPROFRI  ATE  — 0  #  1#  2*-# 

7  '#  OR  3  DIGITS  T(v  jHG  RIGHT  OF  THc  POi'NT)#'/'  TYPE  CARRIAGE  R'# 

S  'ET'JRN  AFTER  TYPING  PAC  H  INPUT  VALUE*  THE  MAXIMUM'/'  NUMBER  A'# 
?  'LL  GW  ED  FOR  EACH  OF  (D#  (2>#  AN  n  (3)  IS  "S"*'/) 

TYPE  199® 

1959  FORMATC  '  0  '  »F  X  #  '  OUTPL‘T  COMMENTS:  THE  ACTION  CRRER  IS  A  SF.OUENCB  !# 
l'OF  NUliDcPS  TO  BF'/«  INTERPRETED  AS  FOLLOWS:  THE  NUH«FR  IN  THF  '# 
2'LEFTMOST  POSITION  OF  THE'/*  SEQUENCE  RE°RESFNTS  that  ACTION  '# 

3 ' TO  9E  7 *KEN  IF  THE  FIRST  FORECAST'/'  CATEGORY  IS  FORECAST*  '# 
5'THE  NUMBFR  IN  THr  POSITION  IMMEDIATELY  TO  THF'/'  RIGHT  OF  '# 

5 'THE  FIRST  POSITION  IN  THF  SEQUENCE  REPRESENTS  THAT  ACTION'/ 

6'  TO  «E  TAKEN  IF  TRF  SECOND  FORECAST  CATEGORY  IS  FORECAST;  «# 

7 ' AND  SO'/'  ON#  UNTIL  THE  NUMBFR  IN  THE  RIGHTMOST  POSITION  OF  '# 

£ ' THE  SEQUENCE'/'  REpRFSE  NTS  THAT  ACTION  TO  BE  TAKFN  IF  THE  •# 
5'LAST  FGRFCAST  CATEGORY'/'  IS  FORECAST  0  THE  ACTION  ORDER  '# 

1 '  THUS  BECOMES  A  STRATEGY  (I*E«#  A'/'  SET  OF  DFCISION  RULES)  '# 

2 '  FOR  THE  DEC  IS  ION-MAKERo ') 

TYPE  1003 


o  o  o  o 


i  <U$  A'F.ETAC  >M  IK'LO  5#  F4  ft  FRI.  il-AU3-72  5*18 AM  F/GE  1 : 2 


100?  FORMAT  C»OROT$.  OF  RUCK  It  !'///> 

TO  TYPE  100,4 

1004  cORM\T  (.'ONjU.  OF  ACTIONS  =  •'/) 

ACCEPT  10  05,  MACTM? 

1005  'FORMAT  (ID 
TYPE  1010 

1010  FORMAT  .(  »  ONQ  o  OF  OSS  CATFG'ORirS  ■  •«/) 

ACCEPT  1 015*  NSTATE 

1015  FORMAT  (ID 
TYPE  102  0 

1P?0  FORMAT  (  •  ON O.o  OF  fq$t  CATEGORIES  =  ‘/> 

ACCEPT  1 0?5 1  N  03$ 

10,05  FORMAT  (11) 

TYPE  1040 

d  0U0  FORMAT  COLOSSES  (INCUUDF  DECIMAL  POINT.):1-/) 

HO  1  1=1, NSTATE 
DO  1  J=i,NACTNS 
TYPE  1044,  I,J 

10«4  FORMAT  (  »  0(  f » I  li  »>  Dlt  <  • ')  =  •/) 

ACCEPT  10«5i  L  OSTA  3(  I j  J ) 

1045  FORMAT  (F4C2) 

1  CONTINUE 
TYPE  ioao 

1070  FORMAT  COCIIME  »POSA*  IL  I  TIES  :  •  /  ) 

DO  C  1=1, NSTATE 
TYPE  1024,1 

1074  FORMAT  (  »OO  jI1,«  )  =  »/) 

ACCEPT  1,075 j  APRICR(I) 

iO^S  FORMAT  (F4b2) 

2  CONTINUE: 

TYPE  ldso 

105Q,  FORMAT  (<?  OCONO  ITIONAL  PR  0FA3IL  IT  IF.S  D/ > 

,00  2  I,si,» NSTATE 
DO  2  ;1  =  1  # NO Q S 
TYPE  i'05«i  IiJ 

1 0e4  FORMAT  ('•'0(l|llj»j'ill^)'  =  •/) 

ACCFPT  1 055 »  BXPTA3(I,J) 

1055  FORMAT  (FB#2> 

2  CONTINUE 

*****  LOADING  RE0UIR5C  FOPMAT  LINF  BASED  ON  THE  NllMBFR  .OF  FORFCAST 
xx***  CATEGORIES# 

00  45  NW D=1 , 5 

KFORMT (NWD)  =  IF  ORMT (NWD,N0PS-1) 

45  CONTINUE 
MN=0 

C  *****  COMPUTING  -THE  Nl'MBFR  OF  STRATEGIES** 


o  o  o 
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NSTRA'T-NACTN.S#*NOPS 
-  MPROUTsN'PTRAT 

c 

,0  QUESTIONS  AN’ 0  I N PUT  TO  RETIMING  lc  SORTED  OUTPUT  IS  DFSIRFQi 

t‘  *****  AN"  IF  SO  NOW  MANY  PTRAtEgIFS  AF£  TO  BE  SORTED® 

C 

TYPE  ;i 075 

1075  FORMAT  (»0D0  YOU  WANT-  THE  STRATEGIES  SORTED  9Y  LOSS#  AND  LI 
1  »■  IN  ASCENDING *  /'  OR  OPR?  (  YES  OR  NO)*/) 

o  ACCEPT  lOFOi  SORTS D 

10*0  FORMAT  (A?)  .  . 

IF  (SORTFDoNS® YES)  GO  TO  «H 
Tvog  1035 

1095  FORMAT  COHOU  M A N v  STRATEGIES  WOULD  YOU  LIKE  LISTED?.*/ 

1  '  <  2  -  n  I R  T  T  I NTCGFR--MAX  IS  50)'/) 

AOCE°T  lOROi  N3R0UT 
1  0°C  F  CRH  AT  (12)  • 

C 

C  *****  COMPUTING  STRATFGY  LOSS  VALUES  USING  EAYCDIAN  APPROACH* 

C.  *****  FOR  DETAILED  EXPLANATION  OF  THE  BAYS  SIAN  METHOD  REFER  TO 
C  *****  ''t:LEMENTARY  DECISION  THEORY"  QY  CHEPNOFF  AND  MOPES  (CHAPTERS  1 
■C  *****  AM0'4)»  PU9LI8HFD  PY  W  ILEYj  COPYRIGHT  1955* 

■t  . 

DO  10  NB'R-liNSTRA7 

CALL  ACNORO  (N BR #N OSS #NA  CTNS#  1.08  ACN) 

•  DO  20  NSsiiNSTATE 
DO  30  N0=1#N0PS 
NOOK  WP  SN  0  P3-N9  +i 
NRRACN*IOPACN(NOf'KWC) 

ACNL0S=EXPTAe(NS#N0)*L0STAB(MS»N8RACN) 
STALOP(NS)~STALCS<NS)+ACNLOS 
30  CONTINUE 

STALOS(NP)-STALCS(NS)*APRTOR(NS) 

STRLOS(NBR)aSTRLO?(NOR)  +  STAlOS(N$) 

?0  CONTINUE 

IF- (SORTS Do NE* YES)  GO  TO  21 
*****  BEGINNING  OF  "PUSHDOWN "  SORT  FOR  SORTED  OUTPUT® 

DO  23  NS1  j.NORO  UT 

TF  (STRL0S(N9R)«LF«STRSRT(N))  GO  TO  26 
IF  (3TRSPTCN)oNEoPo)  PO  TC  2? 

STRSRKN)-STRLOS(NDR) 

M I N  S  R  T  (  N  )  “N  3  R 
GO  TO  21 
23  CONTINUE 
GO  TO  21 

?6  DO  27  NN~MER0UT#N+1^-1 
STRSRT(NM)=STRSRt(NN-1 ) 

MINSRT(NN) =H IN  SRT (NN-1 ) 

77  CONTINUE 

STRSRT(N)=3TRL  0S(N5R> 

MIN3RT  (N)-NSR 


CO 


OOO  OOO  OOO  0,0  0  0-0'  o 


«  /--V  *  <•  ■ 
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**»*'*  RESETT  INS  LOSSPF'iE'Y  OFSFRVATroN  CATEGOPY>TO  ZF*Po 

*1  DO  2*5  NS=ijNlSTATF 
STAL  OS (N ?  )rp  o 
?5  CONTINUE 
10  CONTINUE 

*****  TYPING  OUTPUT  HFADI*GSo 
TYPE  1060 

1  OF  C  FORMATC  0 ’  /27y  j  «  YES  ST  PATEGI ES  » //20X  /«  STRATEGY  '  *7X  *  »  AC  TI  ON  *  » 7X  » 
■i  •  LOSS  •  A2i*>  *  NUN  BE  P'1  j6Y>  •  ORDER-*//) 

IF  (SORTFOoEfto YES)  GO  TO  SO 

*****  NON-SORTSO  OUTDUT  ROUTINE 

DO  4  0  1  =  1  jNST° AT 
NFR  =  I 

•CALL  ACNOPO  <A''9R'jN,:0ES> NA CTNSiTOE AC N). 

TYPE  kPORMT*  N’BR*  f  IODACN(NOSKWD)  jNOSKWO=NOBS  *  1  #  - 1  )jS  TRL'OS(  I ) 

00.  CONTINUE 
GO  TO  95 

*****  SORTEO  QUTPU  T  'ROUTINE 

c0  DO  41  I=1jNSR0UT 
NRR  =  MINSRTU> 

CALL  ACNORO  ( N BRjM OPS , NA CTNS 1 1  06 AC N> 

TYPE  KFORMTj  NBRj  ( I  OP  AC N ( NORK WD > t NOGK WD=NOb St  1 1 -3 > t STRS^T  { 1  > 

41  CONTINUE 

«*««*(*(  QUESTIONING  AS  to  WHETHER  USER  WOULD  LIKE  ANOTHFR  TRY  AT  ]  T0 
?5  TYPE  1 06 P 

106  5  FORMAT  (»ONOULD  V0U  LIKE  TO  TRY  AGAINT  (YES  OP  N0>»/> 

ACCEPT  lO^Oi  AGAIN 
1  070  FORM.AT  (A?) 

IF  ( AGAINoNEoYES)  STOP 
C 

c  *****  ZEROING  OUT  ARRAYS  FOR  ANOTHER  TRYo 
C 

DO  64  I=1>N0SS 
I  03 AON  < I ) =0 
*4  CONTINUE 

00  66  1  =  1  jNSTA  TE 
PTALOS ( I )=0o 
APR  I  OR (I >=0o  . 

00  66  J=1 *N0  3S 
L03TAF(I#J)=0o 
cyPTA3(IiJ)=0o 


ooooooooo  oo  ooooooooooci 
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66  CONTINUE 

00.  67  I  =1 ANSTR  AT'. 

STYLOS (I)-Go  ' 

67  CONTINUE 

IF  (SORTFBoiTEoYES)  GO  TO)  63 
n0  63  IsiiNPROUT 
STR8RT  M)s0.‘ 

.  MINSRT  (I  >“Q. 

6 8  C 0NT  IN JE 
00  to  50' 

END 


***»«*«»*»  mx  «»**.*», it  it**?  »«»*«*  »»  ACN0RD  it#*#*##*#*  kk«*«*x*xn**  **«*«»*** 

THIS  SUBROUTINE  COMPUTES  THE  ACTION  ORDER  FOR  F/CH  STRATEGY 
NUN  PER  AS  RFCUlREDo  PARAMETERS  (IN  ORDER  >  ARE:  STRATEGY  NUMPEF.j 
MUMPER  OF  FOPSCA'ST  CATEGORIES*.  NUMBER  OF  ACTIONS  FOR  I N °!l AND 
ACTION- :0 ROE R  ARRA^Yr  F  OR- ,9ll=f-pU-T-#- 


#V*»KK*(»#HM**H**  *#.#  X  It  It  H  X.X  X  X  M  Xfc  DEFI  NITONS  XXXXXXXXXXXX  XXX# XX  MX  XX  XX  XX  XX 

OV I PND  -  DIVIDEND 

I'OP ACM  -  ACTION  ORRFR  ARRAY 

ICUOT  -  QUOTIENT 

J  -  00-L00P  COUNTER  A'NT  INDEX 

NACTNS  -  NUMBER  OF  ACTIONS 

NFD  -  STRATEGY  NUMBER  ^ 

NOBS  -  MUMPER  OF  FORECAST  C ATF  GO  PIES  .  . 


XX#x»xx«»x#x#x###x-;,x*x##KxpftOGRAM  LISTING  xxi>#*x*x#x  ox  ************ 

SU3R0UTINF  ACN'ORD  <NPR*N  OPSjNACTNS  j  I  OP.ACN) 

INTEGFR  DVIDNO 
DIMENSION  iOPACN(NOeS) 
inUOT=NBP-l 
DO  10  J-liNOBS 
nVlDND  =  IOUOT. 

IF  (IQUOToEQoO)  GO  TO  9 
ID'JOTrOV  ITNO/NACTf'S 
9  I03ACN< j)=OV IONC-NACTNS#  IOUOT+1 
10  CONTINUE 
RETURN 
END 

Va 

C  *»‘XXXX*XXXXXXXXXX*XXX**XXXXXXX.XXX*XttXXXXXXXXXXXX**XX#l'XXXXXXXXXXX#XXX 

*  * 


******  °AVFS  ST& ATEGIE3  ****** 


PAGE  1 


ARE  YOU  EXPERIENCED?  (YFS  0°  NO) 
NO 


7 K IS  PROGRAM  MILL  C-OMPUTE  LONG-RUN  MINIMUM  A.VF.PAGF  LOSSES  FOR 
SF  rC  I  FTC  OPERATIONAL  STRATEGIES  ("SAYES  STRATEGIES")  5  AS  ED  ON  OPERA¬ 
TIONAL  LOSS  VALUES  AND  HEATHER  OCCURRENCE  PRQEAOILITIE So  TH  =  IN°UT 
PARAMETERS  ARE  AS'  FOLLOWS S  <i)  THE  N'UMSFR  OF  ACTIONS  AVAILABLE* 

(2  )  THE  NUMBER.  OF  OR  SERVATI ON  CATEGORIES  USED*  (3)  THE  NUMBER  OF 
FORECAST  CATEGORIES  UScO*  (U)  THE  OPERATIONAL  COST  (LOSS)  cOR  - ACH 
COMBINATION  OF. ACTIONS  AND  OBSERVATION  CATEGORIES j  (5)  thE  CLIMUOLOG 
ICAL  PR03ADI LI TY  OF  OCCURRENCE  OF  EACH  OBSERVATION  CATEGORY*  AND 
<£>  THE  CONDITIONAL  PROBABILITY  FOR  EACH  COMBINATION  OF  OBSERVATION 
AND  FORECAST  CATEGORIES  (1«E®*  THE  .PROBABILITY  THAT  A  PARTICULAR  FORE 
CART  HAS  ISSUER  GIVEN  THAT  A  PARTICULAR  OBSERVATION  CATSGRY 
OCCURRED)  o 

THE  IN  PUT  F.CRM  ATS  APE  AS  FOLLOWS:  CD*  <2>*  ANn  (?)  SHOULD  PE 
TYPFD  AS  ONE  D I  Gil  T  INTEGERS?  ( *0  *  (F>*  AND  (6)  'SHOULD  RF  TYPED  AS 
THPFE  DIGIT  REAL  NUMBERS  (INCLUDING  DECIMAL  POINT  INSERTED  AS 

APPROPRIATE - Oj  1*  2*  OR  7  DIGITS  T0  THE  RIGHT  OF  THE  POINT)® 
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